Reproduced  by 

Armed  Services  Technical  Information  Agenc 

DOCUMENT  SERVICE  CENTER 

KNOTT  BUILDING,  DAYTON,  2,  OHIO 

AD- 


On  the  Singularity  at  a Concentrated  Load 
Applied  to  a Curved  Surface 
By 

E.  Sternberg  and  R.  A.  Eubanks 


A Technical  Report  to  the 


Office  of  Naval  Research 
Department  of  the  Navy 
Washington,  D.  C. 


Contract  N7ora—32906 
Project  No.  NS  035-302 


Department  of  Mechanics 
Illinois  Institute  of  Technology 
Chicago,  Illinois 


November  1,  1953 


On  the  Singularity  at  a Concentrated  Load  Applied  to  a Curved  Surface* 


By 

E.  Sternberg  and  R.  A.  Eubanks 
Illinois  Institute  of  Technology 


This  paper  deals  with  the  singularity  at  the  point  of  application  of 
a concentrated  load  acting  perpendicular  to  a curved  boundary  of  an  elastic 
body.  In  a neighborhood  of  the  point  of  application,  the  boundary  is 
assumed  to  be  representable  by  a sufficiently  smooth  arbitrary  surface  of 
revolution,  the  axis  of  which  coincides  with  the  load  axis.  In  the  event 
the  surface  is  locally  analytic,  it  is  shown  that  the  singularity  is 
identical  with  that  appropriate  to  a concentrated  load  applied  normal  to 
a plane  boundary  if  and  only  if  the  curvature  of  the  meridian  of  the  sur- 
face vanishes  at  the  load  point.  The  required  modified  singularity  for 
the  case  of  non-vanishing  curvature  is  determined  in  closed  form  to  the 
extent  where  the  residual  surface  tractions  are  finite  and  continuous. 

Introduction 

The  singularities  encountered  in  concentrated  force  problems  of  the 
theory  of  elasticity  require  special  and  separate  treatment  if  one  is  to 
arrive  at  practically  useful  representations  of  the  solution  to  such  prob- 
lems. Indeed,  in  order  to  assure  results  which  are  amenable  to  a complete 
numerical  evaluation,  it  is  essential  to  determine  in  closed  form  the 
relevant  singularities  at  least  to  the  extent  where  the  residual  problem 
is  characterized  by  finite  and  continuous  surface  tractions.  This  process 

The  results  presented  in  this  paper  were  obtained  in  the  course  of 
an  investigation  conducted  under  Contract  N?onr-32906  with  the  Office  of 
Naval  Research,  Department  of  the  Navy,  Washington,  D.  C. 
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was  carried  out  in  a previous  paper  by  F.  Hosenthal  and  one  of  the  present 
authors  £l]^  in  connection  with  the  problem  of  the  sphere  under  concen- 
trated loads.  In  this  particular  case  it  was  found  that  the  singularity 
at  a concentrated  normal  force  is  not  identical  with  the  well  known 

singularity  corresponding  to  a concentrated  load  applied  perpendicular  to 

* 

a plane  boundary,  and  that  certain  supplementary  singularities  have  to 
be  introduced  to  effect  a reduction  of  the  problem  to  one  obeying  the  fore- 
going regularity  requirements. 

The  investigation  just  cited  suggests  more  general  questions:  Under 

2 

what  circumstances  does  Boussinesq's  solution  £2]  for  the  semi-infinite 
medium  bounded  by  a plane  and  acted  on  by  a concentrated  surface  load, 
supply  the  complete  singularity  appropriate  to  a concentrated  load  applied 
to  a curved  boundary?  What  are  the  supplementary  singularities  needed  in 
the  event  the  Boussinesq  singularity  fails  to  yield  a regular  residual 
problem?  It  is  the  purpose  of  the  present  paper  to  deal  with  these 
questions  on  the  assumption  that  the  load  is  perpendicular  to  the  boundary 
and  that  the  boundary,  in  a neighborhood  of  the  load  point,  is  representable 
by  an  arbitrary,  sufficiently  smooth,  surface  of  revolution,  the  axis  of 
which  coincides  with  the  load  axis. 


The  usual  conditions  imposed  on  the  solution  to  a problem  involving 
concentrated  surface  forces,  are  as  follows:  (a)  it  must  satisfy  the  field 


equations  of  the  theory  of  elasticity  throughout  the  body  under  con- 
sideration;^ (b)  it  must  satisfy  the  boundary  conditions  for  distributed 


abers  in  brackets  refer  to  the  bibliography  at  the  end  of  the  paper. 
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See  also  £3],  p.  191 • 

^It  is  assumed  here  that  concentrated  forces  at  interior  points,  and 
otherwise  singular  body  force  distributions,  are  absent. 
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tractions (c)  it  must  possess  a singularity  at  each  point  of  application 
of  a concentrated  load  such  that  the  resultant  of  the  tractions  on  any 
surface  enclosing  the  load  point,  and  lying  wholly  in  the  region  occupied 
by  the  medium,  tends  to  the  corresponding  prescribed  concentrated  load  in 
the  limit  as  the  surface  is  contracted  to  the  load  point.  As  emphasized 
in  ^l^,  these  three  conditions,  although  necessary,  represent  an  incomplete 
formulation  of  the  problem  and  fail  to  characterize  the  solution  uniquely. 

A unique  characterization  is  reached  by  considering  the  modified  problem 
in  which  each  of  the  concentrated  loads  is  replaced  with  an  arbitrary 
continuous  distribution  of  surface  tractions  over  finite  surface  elements 
surrounding  the  load  points.  The  solution  to  the  concentrated  force  problem 
is  then  uniquely  defined  as  the  limit  of  the  solution  to  the  modified  prob- 
lem  as  the  surface  elements  are  shrunk  to  the  load  points  while  the  re- 
sultants of  the  distributed  tractions  are  made  to  approach  the  given  con- 
centrated forces. 

This  definition,  which  is  analogous  to  Kelvin's  definition  through  a 
limit  process  of  a concentrated  force  at  an  interior  point, ^ is  natural 
both  on  theoretical  and  an  physical  grounds;  its  usefulness  ultimately 
depends  on,  and  is  confirmed  by,  experimental  evidence,  such  as  that  sip- 
plied  by  Frocht  and  Guernsey  QQ*  The  solution  to  a concentrated  force 
problem  which  is  in  accord  with  the  definition  adopted  here,  automatically 
satisfies  the  three  necessary  conditions  listed  earlier.  On  the  other 
hand,  as  demonstrated  in  £lj,  there  exist  pseudo-solutions  which  meet  the 
three  conditions  cited  but  fail  to  agree  with  the  foregoing  limit  criterion. 

‘‘u  only  concentrated  loads  are  present,  the  solution  must  leave  the 
boundary  free  from  surface  tractions. 

^See,  for  example,  £3^*  ert.  130,  p.  183. 


The  validity  of  Boussinesq 'a  solution  for  the  half-space  under  a con- 
centrated load  is  readily  established  by  applying  the  preceding  definition 
of  a concentrated  surface  force,  say,  to  Cerruti's  solution^  for  the  semi- 
infinite body  bounded  by  a plane  and  subjected  to  distributed  surface 
tractions.  The  pseudo-solutions  to  concentrated  force  problems  discussed 
in  p.]  exhibit  a common  property:  their  singularities  at  the  load  points 
are  of  a higher  order  than  that  of  the  Boussinesq  singularity.  On  the  other 
hand,  it  is  clear  from  similarity  considerations  that  the  singularity  at 
a concentrated  load  applied  to  a curved  surface  must  remain  dominated  by 
the  Boussinesq  singularity;  it  cannot  be  stronger.  This  observation  leads 
to  an  additional  necessary  requirement  applicable  to  the  solution  of  con- 
centrated force  problems:  (d)  the  order  of  the  singularity  at  each  load 

-2  7 

point  must  not  exceed  the  order  r of  the  Boussinesq  singularity. 

The  treatment  of  the  general  problem  to  be  considered  presently  will 
be  based  on  conditions  (a),  (b),  (c),  and  (d).  In  the  special  instance 
of  the  sphere  under  radial  concentrated  loads,  the  correctness  of  the  so- 
lution so  obtained  was  verified  analytically  through  a limit  process.  The 
question  as  to  whether  the  four  necessary  conditions  stated  always  assure 
a unique  solution,  which  thus  coincides  with  the  limit  solution  defined 
previously,  is  apparently  not  easy  to  dispose  of  with  complete  generality. 
Until  a general  affirmative  answer  is  established,  solutions  to  concentrated 
force  problems  conforming  to  the  four  conditions  listed  remain  subject  to 
individual  verification  through  appropriate  limit  processes. 


^See  £3],  art.  16U-166. 
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'r  here  is  the  distance  frcm  the  load  point. 


We  now  turn  to  the  determination  of  a sequence  of  singular  harmonic 
functions  which  is  needed  in  the  subsequent  analysis*  It  is  hoped  that 
this  part  of  the  paper  may  prove  useful  beyond  the  present  application. 


A Sequence  of  Harmonic  Functions 

If  (x,  y,  z)  denote  Cartesian  coordinates*  the  spherical  coordinates 
(r,  0*  Y)  are  defined  through  the  mapping* 


x = r sin  0 cos  /,  y = r sin  0 cos  y,  z = r cos  0 
0 * r < 00,  0*0*7?',  o * y * 2 7f. 


<D 


With  the  notation* 


p = cos  0*  p = sin  0* 
1/2  _ 


(2) 


we  have 

5 2 ^ - 

p - (x  ♦ y ) = rp,  2 = rp  (3) 

for  the  relations  between  the  spherical  coordinates  and  the  cylindrical 
coordinates  ( p , y,  z).  Laplace's  equation  = 0*  in  the  axi symmetric 

case*  for  which  H = H(r*  p)>  becomes 


-2 


H ♦ | H.  + 2-  R - % H = 0. 
it  r r r pp  _2  p 

r 

o 

Equation  (U)  admits  the  product  solutions. 


(U) 
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H(r,p)  = £rn  or  r-n_1J  £p,(p)  or  Q^pjJ  (n  = °»  l*  2»  •••)! 


(5) 


where  Pn  and  are  the  Legendre  polynomials  and  the  Legendre  functions 
of  the  second  kind*  respectively* 


e 

Subscripts  attached  to  functions  which  originally  bear  no  subscripts* 
denote  partial  differentiation  with  respect  to  the  argument  indicated* 

^For  a comprehensive  treatment  of  spherical  harmonics*  see  pQ. 


The  interior  spherical  harmonics  r^P^p)  are  regular  throughout  the 
finite  space  whereas  the  exterior  harmonics  r’n"’^Pn(p)  vanish  at  infinity 
and,  with  increasing  n,  possess  progressively  stronger  singularities  at 
the  origin.  If  we  adopt  the  notation, 

H_n(r,p)  = r"npn-l(p)  (n  = 2>  —>*  (6) 


then 


"-n-l  - - 


-i  a 


Si  H-n 


(7) 


so  that  the  sequence  of  exterior  harmonics  H_n  may  be  generated  through 
successive  differentiations  with  respect  to  a of  the  first  order  exterior 
harmonic  H_^  = r"1. 

For  future  reference  we  recall  that  PR(p)  satisfies  Legendre's 
equation. 


(p8^)  ♦ ♦ i)  pb  - 0 

and  dte  the  recursion  formulas, 

p-n-l  * pn 

(2n  ♦ 1)  pPn  » (n  + 1)  Pn+1  + nPn_1 

P2pn  * 1 - PPPn- 

Furthermore,  we  record  the  special  values 

V0)  = <-1>n  p2n.l(°>  = °> 


(8) 
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> (?) 


v°>  * °» 


^2n*f  l(°)  * <2n+1>  V0) 


(10) 


and 


pnd)  * i,  »;<«  - 


(11) 


U5 ■ — — t — ap_ 

The  arguaent  of  FQ  is  henceforth  assisted  to  be  p and  PQ  * , 


The  Legendre  functions  of  the  second  kind  admit  the  representation* 


<^(P)  = | PnCp)  In  YZp  “ ¥n-l(p) 

in  which 


(n  * 0,  1,  2,  •..)  (12) 


v 2n  “I®  2n *5  « B 2n “ 9 p . 

V 1 = TT  Pn- 1 + 5»Tn-T)  ?n-  3 * 5 . (n-7)  Pn-5  ♦ ••• 


n 


Sip  , P , ¥ . = 0. 

m m - 1 n - a*  -1 


(13) 


Since  (^(p)  has  a logarithmic  singularity  at  p = ± 1,  the  harmonics 
rnQn(p)  and  r”n"'1Qh(p)  in  (f>)  are  singular  along  the  entire  z-axis. 

In  the  problem  under  consideration  we  shall  have  need  for  a sequence 
of  axi symmetric  harmonic  functions  which  is  regular  in  the  half-space 
z i 0 with  the  exception  of  the  origin  r = 0,  where  it  shall  possess 
singularities,  progressively  weaker  than  that  of  = r”1.  Such  a set 
of  potentials  may  be  constructed  through  successive  integrations  with 
respect  to  z of  the  first  order  exterior  harmonic  Ve  thus  seek 

a sequence  of  functions  Hn(r,p)  which  satisfies  the  following  require- 


ments: 


v\  = 0, 

£Hn=nHn-l  <»  - M.-), 

Through  a process  of  induction  we  are  led  to 

Hn(r,p)  =rnTn(p)  (n  = 0,  1,  2,  ...) 


Tn(P}  = Pn<P)  L “ Vn(p), 


where 


L * In  r (1  ♦ p)  ■ In  (r  ♦ z), 


(1U) 

(15)U 


K16) 


(17) 


^ ^The  coefficient  n in  the  first  of  (1$)  is  introduced  for  later 
convenience*  Note  that  (15)  is  analogous  to  (7)* 
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and 


n «■  1 


<*»  - 2<=2n>  Pn  - 2 


<-Un*-  (a.  + U pb 

(n  - mj(n  <f  m ♦ TJ  * 


T.-°> 


(IB) 


provided, 


n 


(n  = 1,  2,  ...). 


(19) 


In  order  to  verify  that  the  functions  Hn,  defined  by  (16)  to  (19) 
have  the  properties  (lU)  and  (1$),  we  first  observe  that  a lengthy  confu- 
tation, based  on  (18),  (19),  and  (9)»  yields  the  recursion  relations, 

P2  vn  + "Py„  ' nVn-l  * Pn> 


n(nVn  - pv;  ♦ »;.!>  - < - 

(n  ss  1,  2 , •••), 

which,  in  turn,  imply  the  differential  equation, 


(20) 


JP  ®2  Yn>  + "(”  ♦ 1)Vn  * 2<pn  " Pn- 
From  (16),  (20),  (21),  (9),  we  have 

* ”PTn  - nTn-l  " - pPn> 
T'(nI„-PTntTn-l)  = 2Pn-l-pV 


(n  * 0,1,2,. ..).  (a) 


(n  = 1,  2,  •••), 


^ (P2T^)  ♦ n(n  + l)Tn  = - (2n  + l)Pn,  (n  = 0,1,2,...)*  (23) 

With  the  aid  of  (16),  (22),  (23),  and  (1*),  Equations  (lit),  (15)  now  readily 
follow.  For  future  purposes  we  list  the  special  values, 

~ ! 9T 

iZNote  that  Tn  * 


9 


8»<°>  ' cnV0)- 

’>>  * P2n(0)-|><0>]  ’X* 

f 

w°>  - [<*  ♦ »V0’] 


-ll 


(2U) 


which  are  obtained  by  induction  with  the  aid  of  (18),  (19),  (22),  (9),  end 

(10). 

According  to  (12),  (13),  (18), 


Vn(p)  - (-l)n  Vn(-p)  = 2In-1(p)  (25) 

and  (12),  (16),  (17),  (25)  yield 

Hn(r,p)  - (-1)“  Hn(r,-p)  * 2rn  Q^p).  (26) 

This  establishes  the  connection  between  the  harmonic  functions  Hq,  which 
are  no  longer  product  solutions  of  Laplace's  equation,  and  the  spherical 
harmonics  of  the  second  kind  given  in  (5)*  Finally,  we  record  explicitly 
the  first  three  members  of  the  aggregate  HR(r,p) : 

Hq  = In  r(l  p),  ^ ■ r Tp  In  r(l  ♦ p)  - ll, 

2 -(27) 

Hj  * \ [(3p2  - 1)  In  r(l  ♦ p)  - p2  - 3P  ♦ l]. 


Sequences  of  Singular  Solutions 

The  sequence  of  harmonic  functions  established  previously  may  be  used 
to  generate  certain  sequences  of  singular  solutions  of  the  field  equations 
of  elasticity  theory.  In  the  absence  of  body  forces,  and  in  case  of 


torsionless  axi symmetry  about  the  s-axLs,  the  general  solution  of  the  dis- 
placement equations  of  equilibrium  may  be  written  as,^ 


20  £u^,  Uy,  uj  ■ grad  (0  + s#0  - £o,0,U(l  - j>)<^J 

(28) 

and 

V^(/>,8)  ® rfy(p»a)  * 0, 

(29) 

13See  C2],  KJ. 
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where 


tv  V ”•] 


are  the  cylindrical  components  of  displacement*  whereas 
Q and  t>  designate  the  shear  modulus  and  Poisson's  ratio*  respectively. 

We  define  two  sequences  of  particular  solutions^  M •»* 
in  terns  of  their  generating  stress  functions  as  follows t 


[Aj  ...  0 = Hn(r,p),  ^=0  (n=  0,1,2,...) 

CBj  ...0  = 0,  V = (n  = 0,1,2,...) 

in  which  Hn  is  defined  by  (16)  to  (19).  The  spherical  components  of  dis- 
placement up,  Ug  and  the  spherical  stress  components  o^*  o^,  oy*  7^g, 
belonging  to  (28),  (29),  were  given  explicitly  in  ^7]].^  Substitution  of 
(30)  into  (b)  to  (7)  of  Reference  Q3*  uae  (22),  (23),  (8),  (9), 
yields  the  following  representations  of  0 W in  spherical  co- 
ordinates . 


r (30) 


20ur  . r"'1(nln  * Pn),  20ifc  = - p r“"  1 
«T  = r”'  1 [n(n  - l)Tn  + (2n  - l)pj 

°e  ■ " 2 [‘  A,  ♦ K - **n] 

°y-  r"‘2  [oTn  - pTn  * Pn] 

7^  a - p rn”  2 [(n  - 1)T^  + pj 


K31) 


^Throughout  this  paper  capital  letters  in  brackets  denote  either  the 
displacement  vector-field  or  the  stress  tensor-field  of  a solution  of  the 
field  equations,  and  equality,  addition,  as  well  as  multiplication  by  a 
scalar,  are  to  be  interpreted  accordingly. 


vanish  identically  in  view  of  the  assumed  rotational 


11 


w« 


20up  * prn  £(n  - 3 + U>>)  TR  + pj 

20iig  = prn  [(3-liV)  Tn  - pT^ 

o'  = r0”1  f"n(n-  3 + 2 >*)pT_  + (2n- 3+ 2 V)pP  - 2^p2t'1 
r l u n nj 

-n(n  + 2^)pTn-  2(n  + >>)pPn+[l-  (3-  2>>)p2]  \ 

(l-2P)p(nTn+Pn)-[^l-(l-2V)p2J 
« - 


— n-1, 

Oq  = r < 


r (32) 


y=  r 

7^  = prn”  1 £(1-  2V)(nTn+  PR)-  (n-  2+  2 ^ )pT^ - pP^j 


Let  Z(r)  be  the  resultant  force  of  the  tractions  to  which  a stress 
field  with  rotational  symmetry  about  the  z-axis  gives  rise  on  a hemisphere 
centered  at  r = 0,  lying  wholly  in  the  region  z ^ 0,  and  having  an 
outer  normal  which  is  directed  toward  the  origin.  If  Z(r)  is  positive 
when  its  sense  is  that  of  the  positive  z-axis,  it  was  shown  in  [8  f6  that 


Z(r)  = 27fr  ■< 


+ 2(1  - >>)r 


(33) 


5r«Jp  = 0 

where  0(r,p)  and  ^(r,p)  are  the  corresponding  generating  stress 
functions.  Applying  (33)  to  (30) > we  find,  with  the  aid  of  (9),  (10), 
(11),  (16),  (17),  (23),  and  (2h),  the  values  of  Z(r)  for  £A J and 

cy- 

27f(-l)n[2n  r213  r n 

[AgJ*  Z(r)  = ^5 — ? [2n  m r + 1 " ^nj 


CA2n  + l3‘  Z(r)=’ 


W 

2*(-l)n  22n(ln)2  r2"4 1 
[2n 


K3U) 


16, 


See  Equations  (31) , (32)  of  Reference  £8]. 


[Bj.  Mr) 


lin'd-  vx-i)”  ■£*  (la)8  r8”*1 


CBJn-J’  Z(r) 


l2n  + 1 

-Urf(l->>)(-l)nl2n+2  r2® 


a8"  (In)4 


[ 


2n  In  r + 1-  2ne 


J 


(31*) 

(Coni.) 


It  follows  from  (3l)  that  the  aggregate  of  solutions  defined  by  the 
linear  combination, 


[Eg  - 2ft  - ><)[intl]  + (n  + l)[Bg 

(n  = 0,  1,  2,  •• • ) 


(35) 


is  self-equilibrated  in  the  sense  that  Z(r)  = 0 for  £eJ.  For  later 
convenience,  we  record  here  the  displacements  and  stresses  of  the  first 
two  members  of  the  sequence  which  also  appear  in  [YJ. 

PJ* 

2Gu  = - (1  - 2*>)pL  ♦ (3  - 2 >>)p  - 2(1  - 7>) 


2GUg  = p T(1  - 2 >*)L  - (3  - 2t>)  Th] 

^.  = -|[(i- 

^ = " 2>)>pJ  - *0 

°y  = F(lVpT [(1  * 8y>p"  8] 

^'hApt[(1-8V)p-1] 


\ (36) 


M*  r 1 

2Gur  = r|2p(l+p2)-p^L+(3->J)p2-2(l4V)p+l-^| 

2Gug  * pr^-2(l-f>9pL+  j^V-3  + >>)p2  - (1-  3>>)p-  2(1-  2V)]J 

ois  -2(l4-^)(p2L+p)4  (5+^)p2-l+ 
r 

Oq  = -2(1  + ^p2L+jJ^-(5+>>)p2-(3->0p+6J 

oy=  -2(i+*>)l+t|-[(-3*i’)p+i-'>’] 

= -2(l+>>)ppL-f^|^  ^-(5+V)p2-  (3  - V)p  ♦ 2] 


(37) 
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Solutions  C*n!]>  [Bj,  £sj  (n  = 0,  1,  2,  • ••),  with  increasing  n, 
possess  progressively  weaker  singularities  along  the  portion  % * 0 of 
the  z-axis,  and  are  otherwise  regular  In  the  finite  space)  for  n * 0 
these  solutions  vanish  as  s -♦  oo  so  far  as  the  stresses  are  concerned. 


Construction  of  the  Singularity  due  to  a Concentrated  Surface  Force  on  a 
Curved  Surface 

Boussinesq’s  solution  £2j^  corresponding  to  a concentrated  load  of 
magnitude  |Ql  applied  at  r = 0 normal  to  the  plane  boundary  of  a medium 
occupying  the  half-space  z ^ 0,  appears  as 

P J - - {»  - * ' * ) C*J  ♦ DLj}  . <38)18 

where  [[A J is  defined  by  (30)  and 


[B_J  ...0=0,  fjJ  = H_^  = r_1. 


(39) 


The  displacements  and  stresses  of  [[Sj  follow. 


”8  = Sm  - 3 ♦ k »*] 


- <X1  - 2 »HpS  ♦ p - I) 

’ ♦ p) 


°b 


2rfr2(l  + p) 
2f^rZ(l  ♦ p) 


KUO) 


Evidently,  2(r)  * Q for  [Sj,  as  is  readily  verified  by  means  of  (38), (39), 
(30),  and  (33). 


^See  also  £3j»  p.  191. 
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If  Q > 0 the  force  has  the  sense  of  the  positive  z-axis. 


Ve  now  turn  to  a concentrated  load  acting  perpendicular  to  a curved 
boundary  (Figure  1).  Let  the  point  of  application  of  the  load  of  magnitude 
|qJ  be  the  origin  and  again  let  Q > 0 if  the  load  acts  in  the  positive 
z-direction.  Suppose  that  the  boundary,  in  a neighborhood  of  the  origin, 
coincides  with  a surface  of  revolution  2I»  the  semi-meridian  P of 
which  admits  the  representation. 


Pi  z as  f (p)  (0  4 p < £) 

f (0)  * f‘(0)  = 0, 


' 

(hi) 
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and  let  f (p)  be  at  least  four  times  continuously  differentiable.  Since 
f 1 (0)  = 0,  5T  possesses  a uniquely  defined  tangent  plane  at  its  vertex 
r = 0.  The  function  f (o)  has  a Taylor  expansion  of  the  form. 


f (/>)  = a 2p2  + a^p3  + 


(12) 


so  that 


a*  » w f"(0)  = 


•2  - ^ = -j-, 

where  k (j>)  is  the  curvature  of 


i3  s|f«»(0) 

P and 


ko 

T 


(1*3) 


ko  ■ k<°>>  K ■ [f J^-  <“*> 

With  a view  toward  examining  the  nature  of  the  singularity  at  the  load 
point  r a 0,  we  first  determine  the  normal  and  shearing  tractions  O'  and 
(Figure  1)  induced  by  the  Boussinesq  solution  W on  the  arc  P . 

To  this  end  we  observe  by  means  of  the  law  of  stress  transformation  and 
elementary  geometric  considerations  that  along  P, 
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Note  that  the  regularity  restrictions  on  f iiqply  the  same  degree 
of  smoothness  for  £ only  if  the  (right-hand)  derivatives  of  odd  order 
of  f (p)  at  j>  - 0 vanish.  Otherwise,  the  even  extension  of  f (j>)  ex- 
hibits discontinuities  in  its  derivatives  of  odd  order  at  the  origin. 

^The  subsequent  analysis  is  not  valid  if  r * 0 is  a conical  point 
of  2* 
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& = 


in  which 


& 4 & - di 

-■  4 — ? -■  cos  2(J3  4 0)  - fa  sin  2<p  4 0) 

o'  — oi 

■r  y £ sin  2(p  +fl)+  ^ cos  2(p  4 Q) 

j3  is  the  inclination  of  the  tangent  of  P . Evidently* 


Kh5) 


cos 


r i + (f  r 


sin 


2f. 


2D=  ^t-h. 

^ 1 4 (t  )2 


m 


Recalling  that 


ire  reach  by  means  of  (3),  (ill),  (1*2),  the  expansions, 

p = a?p  + a^p2  4 (a^  - *^)  p*  ♦ ••• 

*2  2 3 

P = 1 — 5- />  - ^ay5  ♦ •••» 
valid  on  P.  From  (1*6),  (1*8),  and  (2)  follows 

cos  2(p  4 9)  = - 1 4 2ajjjp2  4 8a2a^  4 ... 

sin  2 (p  4 0)  = - 2a^>  - Ua^p2  - 6(a^  - a^)^  4 .... 


(U7) 


KW) 


Kl*9) 


Substitution  of  (1*9)  and  (1*0)  into  (1*5),  after  a lengthy  routine  compu- 


tation, and  final  use  of  (1*3),  yields 

PJ1 


«tf)  » 


Q(1  - 2V)  k‘ 


■ a[ 


TR 

-U(l  - 2>>)k 


4 0(1)  a 0(1) 


4 (5-l»)>)k2-2(l-2>Ok^]  ♦ 0(1) 


(50) 


where  o(l)  and  0(1),  respectively,  denote  functions  of  p which  tend 
to  zero  or  to  a finite  limit  as  p approaches  zero. 

If  designates  the  complete  solution  to  the  concentrated  force 

problem  under  consideration,  we  may  write, 
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1*1  - Eso3  ♦ l*J>  <*> 

in  which  stands  for  the  solution  to  a "residual  problem",  the  charac- 

terization of  which  is  implicit  in  (51).  We  shall  suppose,  without  loss 
in  generality,  that  the  boundary  of  the  medium  is  free  from  distributed 
tractions  on  Furthermore,  we  shall  agree  to  call  J^Rj  "regular" 

on  if  the  residual  problem  on  J is  governed  by  finite  and  con- 
tinuous surface  tractions.  This  requires  for  w* 

<*/>)  = 0(1),  /(/>)  = 0(1).  (52)a 

Equations  (50),  (51)  now  permit  the  following  conclusions.  W *• 
regular  on  ^ if  and  only  if  the  curvature  k and  the  rate  of  change 
of  curvature  k'  of  the  seminaeridian  f of  ^ both  vanish  at  the 
load  point.  We  note  parenthetically  that  k^  is  zero  automatically  if 
the  even  extension  of  P has  at  least  a continuous  rate  of  change  of 
curvature.  In  the  event  that  the  even  extension  of  f (p)  is  an  analytic 
function,  so  that  in  an  analytic  surface,  the  surface  tractions  of 
£rJ  are  analytic  on  if  and  only  if  kQ  = 0.  In  this  case,  and  in 
this  case  only,  the  BoussineSq  singularity  constitutes  the  entire  singu- 
larity of  [3s!]  the  load  point  in  question. 

In  order  to  effect  a reduction  of  the  problem  to  a residual  problem 
which  is  regular  on  V when  k and  k'  are  not  both  zero,  we  need  to 
introduce  supplementary  singular  solutions,  in  addition  to  Moreover, 

in  accordance  with  conditions  (c)  and  (d)  stated  in  the  Introduction,  the 

^Observe  that  /(p)  = 0(1)  is  insufficient  since  it  actaits  a dis- 
continuity in  the  surface  shearing  tractions  of  [[Rj  at  the  vertex  of 

z* 

22 


See  footnote  No.  19* 
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supplementary  singularities  sought  must  bo  self -equilibrated  and  of  a lower 
order  than  r”^.  This  loads  us  to  solutions  W defined  by  (35)  and 
listed  eoqplleltly  for  n * 0,  1 In  (36) , (37)* 

The  local  behavior  of  on  P , In  the  vicinity  of  p * 0, 

Is  established  by  a process  which  Is  strictly  analogous  to  that  followed 
in  the  derivation  of  (50).  We  mayi  therefore,  give  the  results  of  these 


computations  directly. 

PJ* 

°t p)  * -»0  ♦ o(l),  Tty)  ■ j - (1  - 2 i»k0  ♦ o(l) 

w 

dtp)  m 0(1),  y(p)  ■ - 2 ♦ 0(1) 


r (53) 


f (5U) 


It  is  apparent  fraa  (53),  (5U)  that  f of  pJ  still  beoones  Infinite 
at  pm  0 whereas  7*  of  pJ  displays  merely  a finite  discontinuity 
at  the  origin)  the  normal  tractions  & of  both  of  these  solutions  are 
again  finite  and  continuous  on  Z.  as  in  the  case  of  PJ- 

We* now  eet  ^ 


PJ-PJ*PJ 
PJ  ■ PJ  ♦ AJ>J  ♦ \pj, 


KSS) 


J 


and  ask  for  the  values  of  the  coefficients  A > which  assure  the 

0 x 

regularity  in  the  sense  of  (52)  of  the  solution  PJ  to  the  new  residual 


problem.  In  view  of  (50),  (53),  (5U)  this 

^ (l-2>>)Qk 

\ 


to 


|<1+  2(1 -2>»)^] 


(56) 
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supplementary  singularities  sought  must  be  self -equilibrated  and  of  a lover 
order  than  r"^»  This  leads  us  to  solutions  PJ  defined  by  (35)  and 
listed  esqplioitly  for  n * 0,  1 in  (36),  (37)* 

The  local  behavior  of  p J,  pJ  on  P , in  the  vicinity  of  p * 0, 
is  established  by  a process  whieh  is  strictly  analogous  to  that  followed 
in  the  derivation  of  (50)  • Vo  nay,  therefore,  give  the  results  of  these 
computations  directly. 

ot f ) - -2k0  ♦ o(i),  l\p)  ■ i - (1  - 2 P)ke  4 o(l) 

<Kp)  m o(l),  Yip)  ■ - 2 4 0(1) 


It  is  apparent  frcn  (53),  (51t)  that  f of  pJ  still  beoones  infinite 
at  p m 0 whereas  7*  of  pJ  displays  merely  a finite  discontinuity 
at  the  origin)  the  normal  tractions  e*  of  both  of  these  solutions  are 
again  finite  and  continuous  on  2.  as  in  the  case  of  PJ- 
We  ♦now  set 

-PJ*PJ 

pj  - pj  ♦ \pj  ♦ \PJ* 

and  ask  for  the  values  of  the  coefficients  which  assure  the 

regularity  in  the  sense  of  (52)  of  the  solution  PJ  to  the  new  residual 
problem*  In  view  of  (50),  (53),  (5b)  this  leads  to 


(56) 
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WLth  pj  determined  as  in  (55),  (56),  the  surface  tractions  of 
are  finite  and  continuous  on  Jj.  Alternatively,  £s]  admits  the  repre- 
sentation, 

M - PJ ♦ C*J 
W ” PJ  + \ C*o3» 

where  /\q  again  has  the  value  appearing  in  (56)*  The  surface  tractions 
°f  £rJ  are  finite  but  exhibit  a discontinuity  in  f*  at  the  origin. 

It  is  inportant  to  emphasize  that  the  surface  tractions  of  on 

although  finite  and  continuous,  are  by  no  means  analytic  even  if  J 
is  analytic.  This  is  due  to  the  fact  that  the  derivatives  with  respect 
to  p of  o’  and  7*  in  unlike  the  corresponding  derivatives  be- 

longing  to  [S  T,  [Ej,  possess  a logarithmic  singularity  at  j>  = 0. 
Consequently,  the  singularity  of  py,  regardless  of  the  analyticity 
of  does  not  represent  the  entire  singularity  inherent  in  [S]  at 
the  load  point. 

The  sequence  of  singular  solutions  £sj,  thus  G*ves  rise 

to  a sequence  of  solutions  cr<j»  ihi>  cy  to  the  complementary  residual 

problems,  which  are  characterized  by  progressively  increasing  regularity 

in  the  corresponding  boundary  conditions.  The  process  of  successive 

"regularization"  of  the  residual  problem  can  be  continued  in  an  obvious 

manner  by  recourse  to  the  remaining  members  of  the  sequence  of  self- 

equilibrated  singular  solutions  W defined  in  (35) . 

If  ^ is  a portion  of  a sphere  of  radius  rQ,  lying  in  the  half- 

space  z b o,  then  k = r”1,  k'  = 0,  and  the  results  obtained  here  reduce 
r 9 0 0 0 

to  those  given  in  £l]  in  connection  with  the  particular  problem  of  the 

^^The  functions  designated  by  o(l)  in  (5U)  contain  a term, 
g(p)  ■ p In  r(l  4 p). 


4 
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sphere  under  concentrated  loads*  Whereas  the  general  questions  raised  in 
the  Introduction  have  been  answered  within  the  limitations  placed  on  the 
relative  orientation  of  the  load  and  on  the  local  character  of  the  boundary, 
it  should  be  pointed  out  that  certain  difficulties  in  the  large  may  arise 
in  the  application  of  the  results  to  specific  concentrated  force  problems. 

V 

Thus,  the  representations  (51),  (55),  (57)  lead  to  serious  complications 
if  the  line  of  action  of  a concentrated  load  intersects  the  boundary  of 
the  medium  on  opposite  sides  of  the  tangent  plane  at  the  point  of  appli- 
cation of  the  load.  IWe  QjJ,  each  of  which  is  singular 

along  the  negative  z-axis,  give  rise  to  line  singularities  in  the  interior 
of  the  body. ^ finally,  if  the  region  occupied  by  the  medium  is  not  bounded, 
the  representation  (55)  becomes  inadmissible  in  view  of  the  behavior  of 
W at  infinity. 
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